Dynamical Scaling and the Finite Capacity Anomaly in 3- Wave Turbulence 
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We present a systematic study of the dynamical scaling process leading to the establishment of 
the Kolmogorov-Zakharov (KZ) spectrum in weak 3-wave turbulence. In the finite capacity case, 
in which the transient spectrum reaches infinite frequency in finite time, the dynamical scaling 
exponent is anomalous in the sense that it cannot be determined from dimensional considerations. 
As a consequence, the transient spectrum preceding the establishment of the steady state is steeper 
than the KZ spectrum. Constant energy flux is actually established from right to left in frequency 
space after the singularity of the transient solution. From arguments based on entropy production, 
a steeper transient spectrum is heuristically plausible. 

PACS numbers: 47.35.-i, 82.20.-w, 94.05. Lk 



Wave turbulence [1] concerns the statistical mechanics 
of systems containing large numbers of dispersive waves 
which interact conservatively. Such wave systems occur 
in a variety of contexts in nature and engineering. Com- 
monly cited examples include gravity waves on the ocean 
surface, waves in Bose-Einstein condensates and magne- 
tohydro dynamic waves in strongly magnetized plasmas. 
For a short review of applications see [2] . Often the waves 
are driven by external forcing which supplies energy (or 
possibly another conserved quantity such as wave-action) 
at a scale which is widely separated from the character- 
istic scale of dissipation. Since interactions among waves 
are conservative, turbulence results. That is to say, the 
physics is dominated by the flux of energy between the 
forcing scale and dissipation scale. This flux is medi- 
ated by the wave interactions. Theoretically, wave tur- 
bulence is more tractable than its hydrodynamic cousin. 
In the limit of weak nonlinearity, the theory is asymp- 
totically closed under relatively mild assumptions on the 
initial wave statistics (see [3] for a review). Asymptotic 
closure, resulting from the interplay of weak nonlinear- 
ity and dispersion, allows a consistent derivation of a 
kinetic equation describing the long time asymptotics 
of the frequency-space wave action density, N u (t) [18], 
which in turn, determines the leading order behaviour of 
all the higher order cumulants in both Fourier and phys- 
ical spaces. In this article, we limit our discussion to 
the so-called 3-wave kinetic equation (3WKE) describing 
cases where the dominant nonlinearity is quadratic and 
the dispersion relation admits three-wave resonances. 

The 3WKE is a nonlinear integro-differential equation. 
It generally involves 3 scaling exponents, traditionally 
written as a, and d representing the degree homogene- 
ity of the wave dispersion relation, the degree of homo- 
geneity of the wave-wave interaction coefficient and the 



spatial dimension respectively. For isotropic systems, it 
can be written [4] in a form involving only a single expo- 
nent A = (2(3 — a) /a: 



at 



= S^Nu] + S2[N U ] + S 3 [N U ] 



(1) 



where the collision integrals, Si[N u ], S^A^] and 5i[7V w ], 
are written explicitly in the appendix. S^iV^] describes 
forward transfer of energy whereas S2[N U ] and Si[N u ] 
account for backscatter. The details of the wave-wave 
interactions enter the collision integrals through the wave 
interaction kernel, K(uj\,L02) [19]. The exponent A is 
the degree of homogeneity of this kernel. We shall focus 
particularly on the product kernel: 



K(uJi,U2) = (oji w 2 ) - 



(2) 



'Electronic address: connaughtonc9gmail.com 
t Electronic address: anewell9math.arizona.edu 



A key theoretical insight is the fact that Eq. (1) has a sta- 
tionary solution, the Kolmogorov-Zakharov (KZ) spec- 
trum, — ckzlo^ Xkz , which carries a constant flux of 
energy through w-space. The exponent, xkz = -^^j anc ^ 
the constant, ckz, can be found analytically using an el- 
egant technique known as the Zakharov transformation 
[1] . The KZ spectrum is the analogue of the Kolmogorov 
5/3 spectrum of hydrodynamic turbulence. Considerable 
efforts have been made to realise this spectrum experi- 
mentally. Extensive theoretical studies have completely 
characterised its locality and stability properties. 

On the other hand, relatively little is known about the 
time-dependent solutions of Eq. (1). This is the main 
topic of this article. Such solutions are important since 
they describe the process by which the KZ spectrum 
is established when an initially quiescent wave field is 
forced at large scales. Understanding the transient dy- 
namics of kinetic equations like Eq. (1) has become a 
problem of considerable importance in understanding the 
non-equilibrium dynamics of Bose-Einstein condensation 
[5, 6] as well as in wave forecasting. A basic scaling the- 
ory of the transient dynamics was presented in [7] . Two 
important observations were made in that paper. Firstly, 
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FIG. 1: Numerical solution of Eq. (1) with the A = 2 prod- 
uct kernel, Eq. (2), showing the singular transient dynamics 
preceding the establishment of the K-Z spectrum. Numerics 
were performed using the algorithm described in [4]. 

the KZ spectrum is established by the self-similar prop- 
agation of a front in w-space which moves towards large 
uj. Secondly, in finite capacity cases (A > 1) this front 
propagates to u> — oo in a finite time so that the KZ spec- 
trum is set up by a singular solution of Eq. (1). A typical 
example is shown in Fig. 1. When A < 1 - the infinite ca- 
pacity case - the front necessarily leaves the KZ spectrum 
in its wake. For finite capacity cases, the corresponding 
statement is more heuristic for reasons we will examine 
below. The early numerical simulations in [7] indicated 
that the transient spectrum had the KZ scaling. Subse- 
quent numerical investigations of finite capacity cascades, 
firstly in the context of plasma turbulence [8] and then 
in BEC [6] suggested that the transient scaling exponent, 
while very close to the KZ exponent, is not exactly equal 
to it. This fact has become known as the "finite capac- 
ity anomaly" . Investigations of finite capacity cascades 
in differential equations [9, 10] found that the transient 
spectrum is always slightly steeper (for direct cascades) 
than the KZ spectrum, an observation which was heuris- 
tically justified on the basis of entropy production con- 
siderations [10]. 

In this article we perform a careful numerical study 
of the dynamical scaling properties of the 3WKE and 
demonstrate that the finite capacity anomaly is a generic 
feature of finite capacity cascades. The transient spec- 
trum is indeed steeper, although by a very small amount. 
Such a demonstration has been hitherto absent from the 
literature because of the serious numerical difficulties en- 
countered in solving Eq. (1) over a sufficient range of 
scales to measure exponents with sufficient accuracy to 
demonstrate the anomaly. We procede as follows. We 
first introduce the idea of dynamical scaling, define the 
dynamical scaling exponent, a, and show how it relates 
to the transient spectrum. We then show that in the 
infinite capacity case, a is given by the KZ value and 



identify the failure in the corresponding argument in the 
finite capacity case. We then turn to the delicate is- 
sue of numerical measurement of dynamical scaling ex- 
ponents which demonstrates that the transient spectrum 
is steeper than the KZ spectrum in this case. We explore 
how the anomalous exponents depend on the structure 
of the turbulence by varying the amount of backscat- 
ter and the degree of scale-locality of the wave interac- 
tions. Finally we show that a steeper transient spectrum 
is plausible on the basis of a heuristic entropy production 
argument. 

N u (t) tends to a scaling (self-similar) form if there ex- 
ists a monotonically increasing typical frequency, s(t), a 
function, F{x) and a dynamical scaling exponent, a, such 
that 

N^t) ~ S (tf F (^j. (3) 

Here ~ denotes the scaling limit: s(t) — > oo and u> — > oo 
with x = ui/s(t) fixed. The scaling function, F{x), typi- 
cally decays exponentially for large x producing the front 
structure evident in Fig. 1. The small x behaviour of 
F{x) determines the transient spectrum. The general 
properties of the system are well characterised once the 
exponent a and the small x behaviour of F(x) are known. 
We do not know, a-priori, that F(x) is behaves as a 
power near zero. Indeed we already know of one exam- 
ple of decaying wave turbulence in which F(x) diverges 
as x^ 1 log(l/x) [11]. Nevertheless, if we assume power 
law behaviour, F{x) ~ x~ y as x — > 0, and further re- 
quire that the spectrum should remain finite at the low 
frequency end as s(t) — > oo, a reasonable assumption in 
the forced case, then we are lead to conclude that y = —a 
in order to cancel the time dependence. Therefore, all we 
are required to determine is a single exponent, a. 

Substitution of Eq. (3) into Eq. (1) shows that the typ- 
ical frequency, s(t), must evolve according to the equa- 
tion: 

ds 

— = s ? with £ = A + a + 2 (4) 

at 

while F(x) must satisfy the integro-differential equation: 
dF 

aF(x) + x— = StiFix)} + S 2 [F(x)} + S 3 [F(x)}. (5) 
dx 

Three distinct behaviours for s(t) arise depending on 
the value of the exponent £: s(t) grows algebraically 
in time (if £ < 1), s(t) grows exponentially in time (if 
£ = 1) or s(t) exhibits a finite time singularity of the form 

(i* — t)~ f=r (if £ > 1). To determine what actually hap- 
pens, we need to determine a. For forced turbulence, the 
total energy grows linearly in time: J„ u> N u (t) du> — J t. 
Substituting the scaling form, Eq. (3), into this equation 
and differentiating with respect to time yields 

^£ = i -l-a (6) 

dt (a + 2) f °° xF(x)dx ' [ ' 
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Comparing with Eq. (4) fixes a = — ^y^, the KZ value. 
At first glance, we have shown that the dynamical scal- 
ing exponent always takes its KZ value. Care is re- 
quired however: Eq. (6) only holds provided the inte- 
gral J x F(x) dx does not diverge at its lower limit on 

the predicted small x behaviour, F(x) ~ x i~ . This is 
only assured for A < 1 - the infinite capacity case. In the 
finite capacity case, conservation of energy does not con- 
strain the scaling function since the first moment of F(x) 
diverges. Physically, this is not mysterious: the scaling 
solution Eq. (3) does not probe the energy-containing 
scales in the finite capacity case. As a result, the finite 
capacity case exhibits what Barenblatt refers to as self- 
similarity of the second kind [12]: a should be determined 
by solving Eq. (5). 
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FIG. 2: Collapse of the | order moment of the data in Fig. 1 
according to the scaling hypothesis, Eq. (3). The dynamical 
scaling exponent is 2.578 ± 0.021. 

There is no a-priori reason why a should not be given 
by its KZ value in the finite capacity case. Indeed one 
might naively hope it might be based on dimensional 
analysis of Eq. (6). There is, however an established 
precedent to indicate that it probably is not. The kinet- 
ics of irreversible cluster-cluster aggregation (see [13] for 
a review), described by the Smoluchowski kinetic equa- 
tion, has many structural similarities to 3-wave turbu- 
lence, including KZ spectra [14] although it does not have 
an analogue of the thermodynamic spectrum. In fact, 
the Smoluchowski equation simply corresponds to Eq. (1) 
with the backscatter terms, S^A^) and S3(N U ) removed. 
The non-triviality of dynamical scaling exponents is well 
known [15] in that field. An extensive numerical study 
performed by Lee [16] convincingly demonstrated the ex- 
istence of the finite capacity anomaly for the Smolu- 
chowski equation and showed, intriguingly, that it has 
the opposite sign to the putative anomaly in the 3WKE: 
the transient spectrum in the Smoluchowski equation is 
shallower than the KZ spectrum. Our method extends 
Lee's approach to include the back-scatter terms so we 



obviously reproduce Lee's results and various known ex- 
act solutions of the Smoluchowski equation when we turn 
the back-scatter terms off. 

We now turn to the direct measurement of the dynam- 
ical exponent from numerical data in order to demon- 
strate the anomaly in the 3WKE. The details of the nu- 
merical algorithm and its validation are described in [4]. 
There are two principal challenges. The first is the de- 
termination of the typical scale, s(t). The second is the 
determination of the value of a which provides the best 
data collapse according to Eq. (3). Let us now address 
these challenges in turn. The typical scale can be deter- 
mined intrinsically by measuring moments of the wave 
spectrum, M a {t) = J °° lo' 7 N u {t) duj. Eq. (3) implies 
that s(t) is given by the ratio of successive moments: 
s(t) = M a+ i{t)/M a {t). We have already learned to be 
wary of divergences of moments in the scaling limit com- 
ing from the behaviour of F(x) near 0. This issue arises 
again in the definition of s(t). We must take a ratio of 
successive sufficiently high order moments in order to en- 
sure the convergence of the necessary integrals. In this 
paper, we take s(t) = M$(t) / 'M2(t) and restrict ourselves 
to < A < 2 which avoids this pitfall. The challenge 
in determining a from the numerical data lies in the fact 
that it is very close to the KZ value so a very robust and 
sensitive data analysis is required to measure it with suf- 
ficient accuracy. For example, the simulation shown in 
Fig. 1 has a = 2.578 ± 0.021 compared to a KZ value of 
2.5. It is our experience that a "best-by-eye" measure- 
ment of the slope of the wake in a log-log plot such as 
Fig. 1 is completely inadequate. We have instead opted 
to exploit Eq. (3) in its entirety and try to find the value 
of a which best collapses all the curves in Fig. I onto a 
single curve. This was done, as suggested in [17] by defin- 
ing a function which measures the total average spread 
of the collapsed data for a given value of a and then per- 
forming a one-dimensional minimization of this function 
over a. To ensure maximum sensitivity of the analysis, 
the data collapse was actually performed on the logar- 
tithm of (A + 3)/2-moment of N,j(t). In addition from 
removing the subjectivity from the data collapse process, 
the accuracy of the resulting exponent can be estimated 
from the width of the minimum of the data spread func- 
tion. We chose to measure the width of the minimum 
at the 1% level set of the data spread function. A rep- 
resentative sample of the data collapse obtained by this 
method is shown in Fig. 2. It is clearly evident that 
the data collapse is of a very high quality and that F(x) 
is steeper than x -5 / 2 , the KZ value for this particular 
choice of kernel. 

In the numerical simulation, the singularity is regu- 
larised by the onset of dissipation [4] which allows us to 
study the subsequent establishment of the KZ spectrum 
after the transient spectrum has reached the dissipation 
scale. This is illustrated in Fig. 3 which shows the re- 
laxation of the compensated spectrum to the KZ-scaling. 
Note that the KZ spectrum is actually set up from right 
to left. That is to say, the KZ spectrum first emerges at 
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FIG. 3: Relaxation of the transient spectrum to the KZ spec- 
trum post-singularity. 




FIG. 4: Development of a constant energy flux post- 
singularity. 



FIG. 5: Dynamical scaling exponent as a function of A for 
different amounts of backscatter. Lee's results for the Smolu- 
chowski equation [16] are shown for comparison. 




high oj and then propagates to low w as a "backwards" 
front. This interpretation of the post-singularity dynam- 
ics is supported by measurements of the energy flux such 
as those presented in Fig. 4. We remark that the energy 
flux becomes flat as a function of u> first at large cu and 
this region then spreads backwards. The formation of the 
KZ spectrum from the right is consistent with previous 
simulations of the MHD [8] and differential [10] kinetic 
equations although the present numerical scheme is far 
more reliable. This process happens very quickly - notice 
the timings of the successive snapshots in Figs. 3 and 4. 

When it comes to understanding the meaning of the 
anomalous scaling exponents, one may take a mathemat- 
ical point of view and simply view them as being nothing 
more than the exponents defined by the solution of the 
nonlinear eigenvalue problem Eq. (5). From a physical 
perspective this, while clearly a correct explanation, is a 
rather unsatisfactory one and one cannot help exploring 



FIG. 6: Dynamical scaling exponent as a function of A for 
different degrees of locality. 



how the anomaly depends on the physical parameters of 
the problem in the hope of obtaining a deeper insight. 
We performed several systematic studies. Firstly, and 
probably most importantly, we measured the dynamical 
exponent for a range values of A between and 2 taking 
the product kernel, Eq. (2) as a benchmark case. The re- 
sults the circular data points plotted in Fig. 5 along with 
the KZ exponents. They demonstrate that the anomaly 
is generic for the finite capacity regime (A > I) and in- 
creases with A although it is always a small correction to 
the KZ value. From a certain point of view, the small- 
ness of the anomalous correction to the dynamical expo- 
nent is one of its most mysterious features. The other 
sets of data plotted on Fig. 5 show how the anomalous 
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exponents vary with the amount of backscatter. This 
was done by repeating the simulations, first with the Si 
and S2 terms in Eq. (1) reduced by a factor of 2 and 
then with them removed entirely (Smoluchowski equa- 
tion) . While it is clear that the sign of the correction can 
be changed in this way, we stop short of offering a sug- 
gestion of what this actually means. In any case, we are 
not at liberty to tune the structure of the kinetic equa- 
tion in this way in any conceivable experiment. A second 
exploration, presented in Fig. 6, probed how the anoma- 
lous exponent varies with the degree of scale-locality of 
the wave interactions. This was done by introducing a 
deformation of the interaction kernel, parameterised by 
which suppresses interactions between waves of widely 
different frequencies while leaving the overall degree of 
homogeneity, A, unchanged: 



K^(wi,uj 2 ) = K(ui,W2) exp 



-c 



UJ1UJ2 



(u>x + LO2Y 



(7) 

The results show that the anomaly gets bigger the more 
the non-local interactions are suppressed. This was a 
surprise for us but is consistent with the fact that the 
anomalies found in differential approximations [9, 10] 
were somewhat larger. 
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FIG. 7: Sign of the total entropy production on power law 
spectra, Nw = cuo~ x for the product kernel, Eq. (2), for a 
range of values of A. For this kernel, the integral I(x) given 



by Eq. (9) is convergent for x £ 



[4]. 



Finally we offer a heuristic physical argument why 
a steeper spectrum might be natural for the 3WKE. 
Let us formally define an entropy-like quantity, S u (t) — 
log(N u (t)). In the absence of small scale regularization, 
S u (t) diverges on typical realisations of N w (t). Never- 



theless, its rate of production formally makes sense if it 
is defined via the righthand side of Eq. (1): d^*) = 
IvTtT ' ^ ne ma y then ask the question of whether 

the entropy production is positive or negative on a (not- 
necessarily stationary) power-law spectrum, N u = cuj~ x . 
Through application of the Zakharov transformation one 
finds: 



dS„{t) 
dt 



= cuj 2xkz - x - 2 



I(x) 



(8) 



where 
I(x) = 



K( y ,i-y)(y(i -y))- x (i-y x -(i-y) x ) 



(i-y 



2x-\-2 



-(1-2/) 



2x-\-2\ 



(9) 



Whether the entropy production is positive or negative 
for a given spectral exponent, x, depends on the sign of 
I(x). This integral is plotted as a function of x for the 
product kernel with several values of A in Fig. 7. Each 
curve has two zeros, signifying vanishing entropy pro- 
duction (although vanishing for different reasons). They 
correspond to the thermodynamic (x = 1) and KZ ex- 
ponents. The entropy production is always negative in 
between the two. Given that the transient spectrum is 
still exploring its phase space, one would not expect the 
entropy production to be negative. If we therefore limit 
transient spectra to those with positive entropy produc- 
tion, Fig. 7 clearly requires any candidate exponents to be 
greater than the KZ exponent. From this point of view, 
the finite and infinite capacity cases are quite different 
in the sense that for infinite capacity systems no entropy 
production occurs in the wake and the phase space ex- 
ploration occurs entirely at the front whereas in the finite 
capacity case it takes place throughout the full range of 
scales. 

To summarise, we have used what is probably the most 
accurate numerical scheme currently available for solving 
the 3WKE to demonstrate the presence of a small but 
definite finite capacity dynamical scaling anomaly for a 
set of kinetic equations with model wave interaction ker- 
nels given by Eq. (2) and Eq. (7). The results are com- 
pletely consistent with previous observations from the 
original work on the Alfven wave turbulence to the dif- 
ferential equation models suitable for ultra local transfer. 
The natural questions arise: is the anomalous realization 
of the KZ spectrum seen here the general property of all 
finite capacity situations such as, for example, 3D high 
Reynolds number hydrodynamic turbulence? If so, is 
there a general entropy production principle which is re- 
sponsible for guiding the system towards the statistically 
steady state in this anomalous manner? 
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We shall limit ourselves to isotropic, scale-invariant sys- 
tems so that we can exchange the k-space kinetic equa- 
tion for its simpler w-space analogue. 
For the purposes of simplicity, we assume that d = a so 
that there is no difference between the interaction kernels 
of Si [N u ], S 2 [N u ] and Si [N w ]. As a result the thermody- 
namic solution is lo^ 1 . See [4] for details. 



APPENDIX A: COLLISION INTEGRALS AND 
THE WAVE INTERACTION KERNEL 

The explicit forms of the collision integrals which ap- 
pear in Eq. (1) are given below. A full derivation appears 
in [4]. 

Si [N u ] = Jki (u> 3 , uj 2 ) N U2 N U3 5(wi-u)2-u>3) du; 23 

-J K i ( w 3 , OJl ) N u)1 A^ 3 6 (ui 2 - OJ3 - Wi) du 23 {Al) 
-J K i (<*>i , w 2 ) iV Wl N U2 5(uj 3 - lji - u 2 ) div 23 , 



and 



S 2 [N W ] =-[k 2 (lj 3 ,lj 2 ) iV Wl iV W2 5(wi-W2-w 3 )du;23 

+Jk 2 (uj 3 ,lu 1 ) N U2 N U3 5{<jJ 2 -ujz-uji) dw&2 ) 
+ jx 2 {u)i,lo 2 )N Mi N U3 8(u)3-ui-u) 2 ) dcj 23 



S3 [N w ] = -J K 3 (^3 , w 2 ) N Ul N U3 S (lji -uj 2 - 0J3) dcj 23 

+ jK3 {u3,u)i)N Ul N W2 5{u 2 -uj3-ui) dujizki) 

+ jK3 (U!i , 0J 2 ) N U2 N^s S(W3-Ui-U) 2 )(L) 2 3. 



The wave interaction kernels, Ki(uii,u) 2 ) (i = 1,2,3), all 
have degree of homogeneity A = (2/3 — a) /a. There is 
a price to be paid for removing all explicit dependence 
on the spatial dimension from the kinetic equation. The 
kernels appearing in the forward transfer integral and the 
backwards transfer integrals arc not, in general, the same. 
Furthermore, they are not symmetric in their arguments 
as is the case when the kinetic equation is written in its 
usual form (as in [1] or [3] for example) . The relationship 
between them is, however, straightforward: 



K 2 (uii,ujj) = Ki(u>i,i 



K3(u>i,ujj) = Ki(wi,Wj) 



U)i + Ulj 



Uli 



(A4) 



Each of the collision integrals, taken individually has 
a finite flux stationary solution, N u = ckz VJ w ~ , 
where J is the energy flux and ckz is a constant. This 
may be demonstrated by applying the appropriate Za- 
kharov transformations [1] to each integral in turn. The 
stationary thermodynamic solution is hidden in this rep- 
resentation. It is recovered by recombining all three col- 
lision integrals and using the relationships between the 
collision kernels. One obtains the thermodynamic solu- 
tion = ctw^ x ^^~ where ct is constant. In the 
present work, for simplicity, we take a = d so that the 
distinctions between the forward and backward interac- 
tions disappear and the thermodynamic solution is sim- 
ply u)- 1 . 



